FIRST AND SECOND ORDER SEMI-STRONG INTERFACE 
INTERACTION IN REACTION-DIFFUSION SYSTEMS 
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Abstract. Spatial scale separation often leads to sharp interfaces that can be fully localized 
pulses or transition layer fronts connecting different states. This paper concerns the asymptotic 
interaction laws of interfaces for reaction-diffusion systems in the so-called semi-strong regime of 
strongly differing diffusion lengths. A model independent framework jointly for pulses and fronts 
is developped and first order semi-strong interaction is introduced as a general interface interaction 
type. It is distinct from the semi-strong interaction studied over the past decade, which is referred 
to as 'second order' here. Both interaction types are analyzed for abstract systems and for pulses as 
well as fronts. First order interaction for pulses is shown to be gradient-like under conditions that are 
fulfilled for a class of equations including the Gray-Scott and Schnakenberg models. The derivation 
of the two kinds of semi-strong interaction provides necessary conditions on the distinct forms of 
reaction-diffusion systems that can support the interaction types. For illustration, these 'standard 
forms' are applied to analyze a class of two-component systems. 
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1. Introduction. This paper concerns reaction-diffusion systems in one space 
dimension with a scale dichotomy in the diffusion lengths. Such systems are of the 
form 

dtU = Dud.,,U + F{U,V;e) 

dtV = e^D,d,.,V + G{U,V-e). ^'■^> 

with U G K", V G R™, and Du,Dv diagonal matrices with positive entries, and 
X e -D C M with bounded or unbounded domain D. 

The parameter < e ^ 1 is asymptotically small and yields the semi-strong 
limit as £ \, 0. The localization of the V^-components as e \, to a point or jump 
discontinuity defines interface locations x = r^, j = 1,...,A^ of pulses or fronts, 
respectively. The shape of the interfaces is resolved on the small spatial scale ^ = 
{x — Tj)/e. On the other hand, the [/-components remain continuous in a; as e \, 0. 
See Figure [TTTI for an illustration. 

We are interested in quasi-stationary interface patterns, where — > as e — )■ 0. 
The purpose of this paper is threefold. First and foremost, it provides a model- 
independent unified framework for semi-strong interaction of pulses and fronts. Sec- 
ond, it introduces first order semi-strong interaction as a general type of interface 
interaction with Vj = rj{et) for This kind of interaction is distinct from the 

common 'semi-strong interaction', which has rj = rj(e^t). See [5]. It is referred to as 
'second order' here. Third, the distinct equations of motion for first and second order 
pulse and front interaction are derived and analyzed. 

It is shown that first order semi-strong pulse interaction is gradient-like under 
conditions that are fuUfiUed for a class of equations including the Gray-Scott and 
Schnakenberg models. To facilitate the analysis, necessary conditions on F and G for 
semi-strong interaction are derived. These 'standard forms' distinguish the orders of 
interaction and the type of interfaces. However, the existence of interface patterns 
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Fig. 1.1. (a,,b): Snapshots of 4-froi^t solutions to i2.1l on the interval [0,2] with Neumann 
boundary conditions for (a) e = 0.05, (b) e = 0.01. (c,d): Snapshots of 3-pulse solutions to 112.411 on 
the interval [0, 2] voith Neumann boundary conditions for (c) e = 0.01, (d) e = 0.005. The thin lines 
shovi the V -component, the bold lines the U -component. For (c,d) V has been scaled to V = eV. 



that interact semi-strongly is assumed. Sufficient conditions and rigorous proofs for 
the interaction laws are beyond the scope of this paper. 

One may view semi-strong interaction as middle ground between weak interac- 
tion and strong interaction. Weak interaction arises when all components localize 
at interfaces, which therefore interact only through tails that are exponentially close 
to homogeneous equilibria. For instance, a 2-pulse in this case has a spatial pro- 
file that resembles a 2-homoclinic orbit, that is, it twice makes an excursion from a 
saddle equilibrium and passes close to the saddle inbetween. The resulting interac- 
tion law structure is universal and yields exponentially slow motion: in with 
n = roughly rj — rj(cxjp{Kjt/s)), where kj is essentially the slowest spatial con- 
vergence rate to the saddle equilibrium. Such solutions are also called 'meta-stable'. 
See [Tl [TTirni[5^l55] and also Strong interaction occurs, when interface distances 
are on the smallest spatial scale, but to the author's knowledge almost nothing is 
known in this case. An exception is that for scalar equations, such as AUen-Cahn, 
the gradient structure (and other ingredients) sometimes allows to study coarsening 
phenomena of domain walls as in |31| . 'Pulse-splitting' and 'pulse-annihilation' are 
strong interaction phenomena that have been observed in simulations in the semi- 
strong regime [Sl [TUl[2DH^[?71[?8ll30j , but there are no rigorous results and essentially 
no theory. 

In semi-strong interaction the localization of only part of the components defines 
the interfaces due to a separation into large and small spatial scales. In semi- 
strong interaction is driven by the [/-components that do not localize and is hence 
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much stronger than the exponentially slow weak interaction. Indeed, while interface 
patterns in the Allen-Cahn equation interact weakly (are meta-stable), it has been 
shown in [T5] (see also ^2.11 below) that a perturbation of the Allen-Cahn equation 
to the form (jl.ip increases the velocity to order e^. In this paper we show that the 
semi-strong interaction laws are not universal, but generally come in two types: first 
and second order semi-strong interaction. Which of the occurs for a given model is not 
immediately clear, but the standard forms allow to a priori infer the scaling regimes in 
which either type of semi-strong interaction can occur. Moreover, the analysis explains 
order velocities in second order interaction as a result of additional symmetry in 
the limit e — > 0. 

On the other hand, we show that in various models the type of interaction is essen- 
tially determined by the amplitude scaling of a 'feed' term: large feed generates first 
order interaction and small feed the slower second order interaction. Hence, in these 
cases second oder interaction is an asymptotic regime within first order interaction. 

1.1. Relation to literature and models. A dichotomy in diffusion lengths 
occurs naturally when U and V of model densities of particles with strongly 

differing mobilities. Examples in which semi-strong interaction occurs include models 
from chemistry such as the Brusselator, the Gray-Scott model, and the Schnakenberg 
model, as well as the phenomenological Gierer-Meinhardt model for sea-shell patterns, 
and a phenomenological gas discharge model. See [5lll5lll9j and the references therein. 
For scalar U and V , which is the predominant case in the literature, semi-strong inter- 
action has been studied in j5]l71 [?n[2^30ll55] ; an example with two-dimensional U and 
scalar V has been considered in [TnidS] ■ Single fronts have been considered in models 
with two-dimensional V and scalar U in [TTIITSIIMII^ . Semi-strong interaction is not 
restricted to reaction diffusion equations, but has also been studied in a nonlinear 
Schrodinger equation coupled to a temperature field f23| . 

In a large number of models some diffusion coefficients are set to zero due to 
the scale separation. Examples are the famous FitzHugh-Nagumo equations for ac- 
tion potentials in nerve axons (see [H]) and the Orcgonator model for the Belousov- 
Zhabotinsky reaction (see [Hj)- In these cases the interesting interfaces are not singular 
at £ = 0: they do not depend on both the large and the small scale of In ad- 

dition, the interfaces in these models usually move with nonzero speed c in the limit 
£ = 0, and solve the ordinary differential equation in x in the comoving space variable 
X — ct when setting the time derivatives to zero. 

Essentially all previously studied interface motion in the semi-strong regime is 
second order in our notation. Partly as a result and partly as a motivation for the 
present study, in [10] the distinction of first and second order interaction has been 
analyzed in detail for the Schnakenberg model. In our language, [30] studied first order 
interaction in the Gray-Scott model. See also [23 • But it has not been recognised as 
such and the connection and relation to second order semi-strong interaction has not 
been made - the present paper fills this gap. 

There is a fair amount of literature concerning stationary pulse and multi-pulse 
existence and stability in this context. See [51l51 fTCl[T7H2^]l^lK2 6 II 2 8 j and the references 
therein. Additional time scale separation is considered for instance in [51[Tni[H]. The 
existence of stationary pulse patterns in a semi-strong limit has been proven rigorously 
for a class of two-component systems in [5] (and for special models in earlier work 
of these authors), and for a three-component model in [5]. The arising equations of 
motion have been rigorously justified in a certain approximative sense and for special 
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models in [7J[TB] , based on the renornialisation approach developed in . 



This paper is organized as follows. The next section introduces the concepts via 
basic examples. Section [3] gives an overview of the main results. In f^lthe standard 
forms are derived abstractly, and applied to a class of two-component models in fJS) 
The equations of motion for the second order semi-strong case are considered in ^ 
and for the first order semi-strong case in ^JT] 

Acknowledgement. This research has been supported in part by NWO cluster 
NDNS+. The author thanks Matthias Wolfrum, Julia Ehrt and Arjen Doelman and 
Michael Ward for helpful discussions. Special thanks go to Peter van Heijster and in 
particular Tasso Kaper for comments on draft versions. 

2. Basic examples. To motivate and illustrate the subsequent more abstract 
analysis we consider the perhaps simplest models that support first and/or second 
order semi-strong interaction. In the following, subindices denote the order in the 
e-expansions of the solutions, for instance, V ^ Vq + eVi + 0{e'^). 

2.1. Second order front interaction. Let us begin with the following two- 
component model that supports semi-strong front interaction. 

dtU =d,,U^U + V , . 

dtV =e^d,,V + Vil-V^)+eU. ^^■'> 

This model is a perturbation of an AUen-Cahn equation and a case of the FitzHugh- 
Nagumo system with unusual parameters and scalings, and also a reduced version of 
the three-component system studied in [TS]. In Figures fOT a.b). |2. II we plot examples 
of relevant solutions. 

Assuming time-independence to leading order in e implies fixed Vq = V* £ 
{0,1,-1} and 

d^^Uo = Uo + V*. 

This equation describes the (leading order) shape of the 'background field' U and has 
explicit solutions in terms of exponentials. We refer to these as large scale solutions. 

In order to resolve jumps of Vq between the choices for V* ^ it is appropriate to 
switch to the small spatial scale ^ — x/e. On this scale we use small letters u^v and 
(|2.ip becomes 

dfV = d^^v + v{l — v^) + su. ^ ' ' 

In this case assuming time-independence to leading order in e implies constant 
uq (for bounded solutions), and 

diiVo = -vo{l-vl), (2.3) 

which has explicit heteroclinic solutions of tanh-form that allow for connections be- 
tween V* = ±1 in either direction. We refer to these as fronts, and they resolve the 
jumps to leading order. Compare Figure [TTlT a.b). 

We also infer that fronts connecting ±1 must be stationary on the ^-scale: a 
travelling wave ansatz ^ — >■ ^ — with velocity c introduces a friction term —cd^v 
in (|2.3p . Since (|2.3p is Hamiltonian and v ~ ±1 have the same energy, heteroclinic 
connections for c 7^ cannot exist. Therefore, fronts move at most with velocity 
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Fig. 2.1. Slagged V -profiles of the evolving 4--front solution from Fiaure lTTlV a. b ) for (a) e = 0.01 
and up to time 2000, (b) e = 0.005 and up to time 8000. The front locations move essentially on 
the same trajectories in both cases, which corroborates the time scale of second order semi-strong 
interaction. At the coarsening event near x = 1.6, where two fronts annihilate, the asymptotic 
description breaks down. 

on the x-scale and thus system (j2.ip can only support second order front interaction. 
See Figure [O 

A sequence of N fronts generates a leading order solution if the large scale so- 
lutions between these match appropriately at the front locations (and the bound- 
ary conditions if present). This yields a system of algebraic equations, analogous 
to that in [TB], whose solutions are (locally) parametrized by the N front locations 
rj,j — l,...,iV. Accordingly, the linearization of (|2.ip in such a solution has, to 
leading order, a kernel of dimension N spanned by the front translations. This in 
turn gives solvability conditions for a leading order construction of a time-dependent 
solution and thereby provides the laws of motion. Notably, here the term ell in the 
V^-equation comes into play. In SJB] we provide details of this procedure. 

2.2. First and second order pulse interaction. Let us now turn to a simple 
example where first order interaction occurs, namely 



This is a simplified Schnakenberg model, cf. SjH a-nd closely related to the 'linear' 
model in [SDj. See also EH]- 

As before, we assume time-independence to leading order in e. This gives Vq € 
{0, 1/C/o} and 



Possible interfaces are again resolved on the small scale ^ = x/e, which, to leading 
order, gives constant uq and 



dtU =d^^U + a-V 
dtV = e^da^xV -V + UV^. 



(2.4) 



dxxUo = -a -I- Vq. 



(2.5) 



Since the latter only allows for homoclinic connections to vq = (explicit cosh- form), it 
follows that (|2.4p only supports pulse interaction. In contrast to the previous example 



(a) (b) 

Fig. 2.2. Slagged V -profiles of an evolving 3-pulse solution to 112.411 for a = 2.95 and the initial 
condition plotted in Fiaure Vl.lY c.d). (a) e = 0.01 over 200 time units, (b) e = 0.005 over 400 time 
units. The pulse locations move essentially on the same trajectories in both cases, which corroborates 
the e time scale of first order semi-strong interaction. Also the pulse amplitudes are similar, which 
reflects the scaling of V . 




Fig. 2.3. Stagged V-profiles as in Figure V2. 2\ for a = 6, £ = 0.005 and up to time 400. At the 
pulse-splitting event near x = 0.3, the asymptotic description breaks down. 



this means Vb = 0, which imphes that the large scale solution Uq is independent of 
Vq and therefore decoupled from the pulse locations. The result is motion of pulses 
that is driven by a fixed external field given by Uq- Let us have a closer look at how 
this degeneracy arises. Matching of small and large scale for derivatives means that 
d(^^u ~ e^(u — a) should be written as a first order system in the form 

T ^ (2-6) 

Now, matching at a pulse location x — requires that the leading order large scale 
derivatives \ims\o dxUo{x^ ± 5) equal the hmiting leading order small scale deriva- 
tives limj_i.-|-oo p(^). From the equation for p in (j2.6p we infer that the assumption 
limj_j.o — causes the decoupling. 



Therefore, let us make the ansatz V = eV in (|2.4p . which, for leading order 



stationary solutions, gives 



e^d^^V =V-e-^UV\ ^ ' 

Substituting expansions of V and U in e and comparing terms of equal order we find 
the solvability conditions Vo = and Vi G {0, l/[/o}- We shall argue below that 
Vi = so that the leading order large scale problem still reads 

a,,f/o = -a. (2.8) 

However, on the small scale p.Sp turns into 

Q » - -1-2 
O^^W = V — £ UV , 

so that the leading order term uq must vanish and we therefore set u = eu (but not 
U — eU). Hence, changing (j2.7p to the small scale gives to leading order 



%z)o = Vo - uoWo- 



(2.9) 



Recall that matching of small and large scale at a; = means Uo{r^) ~ uq = and 
P± lim d^U{n ±6)^ lim d^U{^), 

(5—^0 ^±oo 



which are thus asymptotic boundary conditions for (|2.9p that are completed by 
lim^^ioo "20(0 — for pulse-solutions. Recall that Vi = I/Uq appeared as a pos- 
sible solution above, but is ruled out because C/q = at pulse locations. 

Let us consider (j2.9p as a 4-dimensional first order ODE with po — d^uo and 
Qo ~ d^vo- It has the two-dimensional invariant space {vq = ~ 0} consisting of 
affine uo(C) = ^0(0)^ -I- 'Uo(O). It is straightforward to compute that this space is 
normally hyperbolic and each point has one-dimensional stable and one-dimensional 
unstable manifolds. Pulse-solutions lie in the intersection of its three-dimensional 
(center-) stable and unstable manifolds, which in four dimensions is generically a two- 
dimensional transverse intersection. This means that we can expect at best a curve 
of P-|--values for which pulse-solutions exist. In fact, due to symmetry, this curve is 
{P+ -I- P_ = 0}. However, the slopes of the large scale solutions already exhaust the 
two-dimensional P±-space so that (|2.9p does not provide sufficiently many solutions. 

The problem is that (|2.9I) enforces motion of at least order e on the ^-scale, which 
is appropriate for second order semi-strong interaction. In the present case, however, 
we have to look for first order semi-strong interaction. This is done by allowing for 
a comoving frame ^ — > ^ — ci before reducing to the leading order equations, which 
turns (113 into 

^ ~ ^ ^ ~ ^2 (2.10) 

Indeed, the resulting motion is first order semi-strong as corroborated in Figure l^?^ In 
Figure fTTlT c.d) we plot the solution profiles, which illustrate the large and small scale 
separation discussed here. Figure 12.31 shows the typical pulse-splitting phenomenon 
in equations such as (|2.4p for large 'feed' a. 
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For small a the phenomenology is more akin to the front interaction of the pre- 
vious example and involves much slower pulse motion as well as coarsening by loss 
of pulses. This motivates to look for second order interaction for suitable a = o(l). 
Using the abstract 'standard forms' derived in ^it is straightforward that the correct 
scalings are a = s^^'^a, U = s^^'^U, V — £~^/'^V. See ^ 

Substituting this ansatz turns (j2.8p into 

and (P?TU1) into 

%uo = 

The latter is again Hamiltonian, so that c = is required for homoclinic solutions, 
which implies second order semi-strong interaction as in i j2.1l 

3. Main results. In this section we summarize the main results for (jl.ip . Our 
basic standing set of assumptions is as follows. 

Hypothesis 1. There exist ytaso yfront ^ ^ interface locations 

Ti < . . . < r N in int(I3), and a family of solutions {U^, Ve) to (jl.ip with the following 
properties for e £ [0,eo), x ^ D and j = 1,...,A'^. The family is bounded, twice 
continuously differentiahle in {e,x) ^ {0,rj), and quasi- stationary: dt{Ue,Ve) = 0{e). 
The limiting Uq is continuous and non- constant m D while Vo{x) G {ybase yfrontj 
and dtVe — 0{£'^) for x ^ rj. As e — > the spectrum of the linearisation of (jl.ip in 
this family consists of an N-fold zero eigenvalue and the remaining spectrum is stable 
and uniformly bounded away from the imaginary axis. 

This paper is not concerned with sufficient conditions under which Hypothesis [1] 
holds. See [71IT6] for (implicit) proofs for certain models. 

Here and in the following we write x± for the left and right limits lim^x^o x ±8. 

Definition 1. For solutions from Hypothesis\^ 

1. An interface is called a pulse if Vo(r'j±) = ^basc ^ front if |Vo(?'j+) — 
Vo{rj-)\ = 1^'^^'^'= - yf™nt|. 

2. The interfaces interact k-th order semi-strongly if there is non-trivial R : 

such that ^r^ = e''R{ri, . . . , tat) + o{e''). 

Starting from Hypothesis [TJ the strategy is to identify requirements on the form of 
F and G so that interfaces can interact first or second order semi-strongly. Notably, 
solutions are assumed to remain bounded as e — > 0. This circumvents searching 
for a scaling oi U, V and parameters that achieves boundedness and semi-strong 
interaction, which is the standard ad hoc strategy when given an explicit rational 
form of F and G. 

We refer to l/^aso ^j^g background state and, to ease notation, assume without 
loss of generality (by shifting V) that y^asc _ 

3.1. Standard forms. The necessary e-expansion and algebraic form of F,G 
for semi-strong interaction derived in 21 reads 

dtU = D,A,U + H{U,V;e) + e-^F%U,V) + e-^FHU,V)U)V 

dtV = e^D^d^^^y + eE{U,V-e) + {G'{U,V)+e~^G^{U,V)U)V, ^ ' 

where E, H are smooth at e = and for fronts F^ ,G\F^, G^ have the factor ]/— 1/^'°"*: 

F%U, V) = F^'iU, V)[V - F^™"'], G"'([/, V) = G'^iU, V)[V - y^""*], 

F^U, V) = F^U, V)[V ~ -l/f™"t], ([/, V) = G^iU, V)[V ~ F^™"*]. 
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In addition, for V* e {0, T/f™nt|^ Qf(^jj^ ^ then also E{U, V^*; 0) = 0. 



Nontrivial second order interaction for symmetric pulses requires F'^ ^ 0, and 
otherwise nontrivial interaction is driven by E as in the example (|2.ip . We refer to 
with i^^ = = as the standard form for second order semi-strong interaction. 
Indeed, all models where second order semi-strong interaction has been found have 
this form. Compare the example (|2.4p and ^J5] In particular, setting ^^ = 6*^ = 
covers the 'normal form' for semi-strong pulse interaction in two-component models 
proposed in [5]. This has n = m = 1 and 

U = -Mf/, = 9{U)V F'^ f{U)V, E = 0. 

Under non-degeneracy assumptions on /, g this form is in fact sufficient for existence 
of pulses as in Hypothesis [T] See [5] . Also the three component model studied in [T5] , 
which has fronts, possesses the standard form for second order semi-strong interaction 
with U = (wi, W2) e M^ F e R, yf""* = 2 (for F in [H] take 1/ = ^ - 1) and 

= 0, ^{awi +I3w2+-f), = F - 1 , H = {{V - wi)/t, {V - W2)/0). 

Concerning first order interaction, note that the standard form covers the example 
(|2.4[) with F^ = 0. Compare fJ5] First order interaction for fronts has not been found 
in a specific equation to the authors knowledge, but arises naturally from the analogy 
to second order front interaction. 

3.2. Lav^rs of motion. The arising laws of motion for interface interaction are 
derived in 21 ^-nd fJBl In this section we summarise these results. We immediately note 
that, due to spatial translation symmetry, single interfaces on periodic D or D = 
have constant leading order velocity, and reflection symmetric pulses are stationary 
due to reflection symmetry of (jl.ip in x. 

On the 'large' x-scale, for x G {rj,rj-^-l), we write solutions Ue,Ve from Hy- 
pothesis [U as Ue = U-i = Ua.j + eC/ij + 0(£^), K = Vj = Vqj + eVi.-j + 0{e^) for 
j = 0, . . . , 4- 1 with boundaries at tq, r^r+i, if present. On the 'small' spatial scale 
~ (x — rj)/e we use small letters u, v and omit the index j on ^ in the following. 

3.2.1. Second order semi-strong interaction. In this case the interface at 
X = rj is a heteroclinic or homoclinic orbit wo.j of 

£'d%woj + ^"'(aj, uoj)uoj = 0, (3.2) 

where aj :— uqj is constant in ^ and thus acts as a parameter. The next order in u 
is given by 

- ^i..(0) +Po,j(0)e - ['^ D-^F%a,,vo,jmvo,jiri)dvdC, (3.3) 
Jo Jo 

and for t = e^t the equations of motion read to leading order 

^r.j = -(5„G'(aj,-(;oj;0)[Mi,j,?;oj] + £^(aj, woj; 0), Wj)/(a{t;oj, Wj). (3.4) 

Here Wj span the kernel of the L^-adjoint of 

Lj = Dyd^^ + G'iaj, vo,j) + dyG'iaj,voj)vQj. 
9 



Note that her Cj is spanned by the translation mode d^voj. In practice, if aU interfaces 
are pulses then vqj = wo,i, and if all are fronts, then wo,2j+i = I'o.ij i^o,2j = ^0,2- 

At regular points of the matching problem between x- and ^-scales, the right hand 
side of p.4p yields a vector field and generalizes the equations of motion reported in 
the literature. Together with matching the IZ-components, the leading order problem 
at one interface typically involves all others, which means that interfaces are globally 
coupled to leading order. 

For TO = 1 (|3.2p is Hamiltonian so that heteroclinic solutions whose asymptotic 
states have different energy (as in [TTIIIHIIMIIIS] ) cannot be stationary. However, two 
equilibria having the same energy is not structurally stable, so that locking motion at 
order for fronts requires additional structure in C, such as symmetry in u and v. 
A trivial case is when G^{u,v) is independent of u as in ^ 32. II and [15] . 

Self-adjoint C and symmetry of the interfaces allow to simplify the equations of 
motion and distinguish fronts and pulses. For instance, for even pulses the term in 
E vanishes, while it essentially drives the motion of odd fronts in case of further 
symmetries. See t)6.2l It is one of the strengths of the model independent approach 
that it provides a common framework, linking the results for different models from 
the literature. 

3.2.2. First order semi-strong interaction. In this case the interface at rj 
and its leading order velocity Cj are determined simultaneously by the boundary value 
problem 

lim vo,jiO = (3.5) 

{— >-±oo 

lim = d^Uoirj±), 

{— J-±oo 

where uqj = and this typically vanishes; else it must be a common root of . 
Note that = = 0, Cj = in (j3.5p gives p. 21) . that is, second order interaction. 

An interface is a solution of (|3.5|) that is homoclinic (pulse) or heteroclinic (front) 
to the invariant manifolds at vq G {0, which consist of affine ui. Compare 

the basic example (|2.4p . To leading order, matching involves only the nearest left and 
right neighbors, because uqj is typically one fixed value. It is useful to view p.Sp as 
an a priori description of the velocities Cj = C {P~ , P^' , u^j) with = d^Uoirj±) as 
parameters. In terms of C and for t = et the leading order equations of motion for 
first order interaction are 

DudxxUo = -J(C/o) , X GUj^Q^Nirj,rj+i) 

Uoirj) = 4' J^l,---,N (3.6) 
= C(aeC/o(r,-),a_eC/o(r,+),4), 

where J{Uo) = H{Uo, 0; 0) - Fl{Uo)E{Uo, 0; 0) for certain F^ (see g]), and boundary 
conditions apply at ro,rN^i. 

A benefit of model independent equations of motion is that for specific equations, 
where numerical observation suggest semi-strong interaction, the laws of motion can 
be readily computed and compared with the simulations. Moreover, general properties 
can be identified a priori. Assuming solvability of the first order equations of motion, 
for scalar u and certain properties of J and C, we prove in fJ7] that first order pulse 
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motion possesses various Lyapunov functionals: in particular, the largest amplitudes 
of U and distances between pulses decay. This severely constrains the possible leading 
order dynamics. For instance, periodic interface motion is not possible. 

4. Standard forms for semi-strong interaction. In this section we derive 
(|3.1|) from under Hypothesis [TJ As a first step, we note that non-constant Uq 

and constant Vq away from interfaces implies that G^, in (jS.ip -j must have a factor 
V -V* for V* = V'°'"'° = in case of pulse and V* € {0, ^f™"*} in case of fronts. 

This a priori rules out semi-strong interaction for equations with E — and G 
linear in U, such as the aforementioned FitzHugh-Nagumo model, where G{U, V) = 
p{U — ^V) (see [H]), and the Oregonator model, where U = (wi, W2) and G{U, V) = 
wi~V (see [9]). 

In the following we consider solutions Ue,Ve from Hypothesis [1] and omit the 
subindex e for brevity. Hypothesis [T] implies that, away from interfaces, 

Dud^.U + F{U,V-e) = 0{e) 
e^D,d,^V + G{U,V-e)V + eE{U,V-e) = 0{e^), ^^■^> 

where the right hand side contains the time derivatives. As before, we write U, V 
on the x-scale and u, v on the ^-scale. We repeatedly make use of the smoothness 
assumption in Hypothesis [T] without explicit mentioning. 

The small scale problem blows up an interface region to the real line via the 
scaled variable = [x — rj)/e, for each j = 1, . . . , A^. By Hypothesis [U we can write 
rj = rj{T, e) with r = et. Setting Cj = -^rj (0, 0), this means = x/e — Cjt + 0{e). 
We consider a single interface and drop the subindex j in the following. Substituting 
the small scale ansatz gives. 

Dyd^^v + cd^v + G{u,v;e)v ^ 0{e). ^ ' ' 

Pulses are solutions of these equations in the limit £ = whose z;-components are 
bi-asymptotic to {v = 0}, and fronts heteroclinic connections between {w = 0,w = 

^frontj^ 

In preparation for a subsequent matching of slopes on large and small scales, as 
in (j2.6p we write (14. 2p as a first order system in the form 

d^u = sp 

Dud^p = -e^cp-eF{u,v;e) + 0{e^) 

d^v = q 

Dyq^ = —cq — G{u,v;e)v + 0{e). 

This choice of p makes the x and ^ derivatives of u and U directly comparable. 

Hypothesis [T] yields an expansion u = uq + 0{e), p = po + 0{e), v = vq + 0{e), 
q = qo + 0{e), u = Uo + O{e). Setting e = in (|4.3p gives constant uo,pq, which means 
that the amplitudes and slopes, respectively, of Uq to left and right of an interface 
are equal to leading order. Since (|4.ip i yields a second order ordinary differential 
equation for Uq, it follows that a pulse interface does not affect Uq. It can thus be 
any fixed solution to (|4.ip i at e = that satisfies the boundary conditions. Hence, 
to leading order, pulses would move in a fixed external field given by the large scale 
solution. However, at a front interface the value of V in (|4.ip jumps between V^^^'^ 
and yfi^o"' and it thus effects the large scale solution as in !j2.1l 
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As in H2.21 the decoupling for pulses disappears when F is of the form F{u, v; e) = 
F{u,v; e). It turns out that the exact nature of the singularities of F and G is 
essential for the type of semi-strong interaction. Recall that here we assume bound- 
edness of U, V and arrive at the conclusion that the right hand side of (|l.ip must have 
a singular term. In contrast, in ^2.21 we scaled v a posteriori to obtain boundedness 
and thereby introduced a singularity in the originally regular right hand side. 

4.1. Second order semi-strong standard form. Let us assume now the sim- 
plest singularity, namely 

F{U, V; e) = H{U, V; e) + e-^F^{U, V), 

and regular G so that G(C/, V: e) = G^{U, V) without loss of generality, by modifying 
E, if required. 

In order to solve on the large scale, we need to evaluate at £ = 0. Since 
V = Ois) for X ^ rj (sec Hypothesis [T]) , this is possible if Fi{U,V) = F''{U,V)V. 
The leading order large scale problem thus reads 

Dud^^Uo + H{Uo, Vo;0) - F,^(C/o)S(C/o, l^o; 0) = 0, (4.4) 

where F^ stems from substituting the leading order correction to the root of order e 
of the ^^-equation G"(C/, V)V + eE{U, V; e) = 0{e^) into F^ Note that here the as- 
sumption dtV = 0{e^) away from interfaces is used (see Hypothesis [T]), and solvability 
requires E{Uo,0;Q) = whenever G^{Uo,0) = 0. The correction F^ in fact vanishes 
in all cases treated in the literature, where E = or F^ = 0, or F^ is quadratic in V. 

The leading order small scale problem (j4.3|) now has the form 

d^UQ = 

Dud^Po = -F''{uo,vo)vo 

d^^vo = qo 

Dvd^qo = -cqo - G''{uq,vo)vo. 

In order to obtain a complete leading order solution, the large and small scale 
solutions need to match appropriately at the interface for x = rj. Hypothesis [1] 
implies Vq = V^^° = on both sides of the interface for pulses. Since uq is constant, 
continuity of Uq and Vq requires uo(C) = Uo{rj) and vo(C) — > as ^ — > ±oo. For the 
derivatives of U and u, matching means to leading order (note vq ~ voj,po — Po.j) 

dMrj±)^Po{0)T F'{Uo{rj),vo)vod^. (4.6) 

Jo 

Interfaces thus correspond to solutions of (|4.5p that are homoclinic to V^'^^'^ = in 
the (wq, 9o)-cquations of (j4.5p for pulses. Similar to the example discussed in H2.21 for 
c = the (vo, (7o)-equations are reversibly symmetric (by reflection symmetry x — >■ —x 
of (jl.ip ) also when including the matching conditions. Hence, homoclinic orbits that 
are reflection symmetric about the wo-subspace are generically robust (codimension 
zero), and persistent under perturbations of mo(0) = U{rj). See [3]. Thus, c = 0{e) 
so that motion is order (at least) £^ on the x-scale. See fjG] For m = 1 the equations 
arc in addition Hamiltonian so that c = 0{e) for any homoclinic solution. 

The situation for fronts is essentially analogous. See tj6.2l for the effect of symme- 
tries in the laws of motion. 



12 



4.2. First order semi-strong standard form. Let us now consider the case 
of other singularities in F, G. From the above discussion, first order semi-strong 
interaction arises if the reversible symmetry is broken so that c = C(l) is typically 
required to locate solutions. For clarity of the exposition, we now focus on pulses. For 
fronts all requirements a.t V = V^'^^'^ equally apply to = yh-ont ^^^^ explicitly 
noted in ^ 

4.2.1. Small scale problem. The next order for singularities is = 0{e~'^) or 
G = 0{e~^). We start out by considering F and suppose there is bounded F{U, V; 0) 
such that 

F{U,V;e)^e-^FiU,V;e). 

Expanding F in e and arguing as in the second order case this means that 

F{U, V; e) = F2{U, V)V + eF^U, V)V + e^H{U, V; e). 

From (j4.3D i it follows that 9^wo = so that uq = is constant. Now the expansion 
w = ?io + sui + ©(e^) in (|431) gives 

Dud^p 
Dyd^q 

Since vq is non-constant in the interface, to compensate the singular coefficient of F2 
requires a root at uq so that for constant uq = we have F2{u, v) = F^{u, v){u — u^). 
This implies that at any interface location the value of uq must be at a root of F^ . 
Indeed, in the example (|2.9p the unique root of in w lies at zero. Without loss of 
generality, by shifting u, we can assume = at least for one interface. 

Setting e = in the equation for q generates a right hand side that is independent 
of Ml and Po, so that the symmetry argument from the second order case applies, which 
means interaction with motion of order - at least for symmetric pulses. Hence, 
generally we need to allow for a term of order in G, and, as for F, this must have 
a factor [u — u^). In particular, must be a simultaneous root of F^ and G^ . Indeed, 
in the example (|2.9p the unique root of G (and F) in u lies at zero. Since terms of 
order e can be absorbed into E, we obtain 

G{u, V- e) ^ G'{u, v) + e-^G^{u, v){u - u^). 

Substitution into (|4.7p gives to leading order 

Dud^Po 
d^vo 
Dyd^qo 

This is the generalisation of (|2.9p and analogously {vq = O^vq = 0} is an invariant 
subspace which consists of affine mi(C) = ei^ -I- 62 for any ei, 62 <E M". This space is 
also the center manifold of any equilibrium with vq = 0, and for c = its transverse 
eigenvalues are the square roots of those of the matrices G^{u^ , 0) + G^{ui, 0)[ui, •]. 
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= —F^{u^ + £Ui,v)v — e ^F2{u^ + eui,v)v + 0{e) 
= q 

= —cq — G{u^ + eui,v;e)v + 0{e). 



(4.7) 



= Po 

= -F''{u^,vo)vo- F^{u^,vo)[ui,vo] 

= go 

= -ego - G^'iu^, vo)vo - G^{u^,vo)[ui,vo]- 



(4.8) 



Matching large and small scale requires po{ioo) = dxUo{rj±)^ but in contrast 
to the second order small scale problem, here 9^po depends on ui to leading order. 
Compare t j2.2l In particular, asymmetry of dxU^ at Vj implies asymmetric boundary 
conditions pq{—oo) ^ ~pq{oo), which break the reversible symmetry at c = 0. There- 
fore, typically c 7^ is required to locate a solution that is homoclinic or heteroclinic 
in the WQ-component. 

4.2.2. Large scale problem. Regularity of (j4.ip i with 

F{U, V) = e-'^F^U, V)V + ©(e^^) 

at £ = generally requires Vq — Vi = 0. However, many examples have F^{U,V) = 
F^ {U,V)V'^ . See SJS] In this case Vo = is sufficient for regularity. We therefore 
distinguish F^{U,0) 7^ from F^{U,0) = and refer to these as the 'linear case' and 
the 'quadratic case', respectively. 
Recall (|IJ|) 9 reads 

e'^E{U, V; e) + eG^iU, V)V + G\U, V)V = 0, (4.9) 

so that Vb = is indeed always a solution. 

The linear case. This case requires Vi = 0, which means for (|4.9p that G^([/o, 0) ^ 
or else E{Uq,Q;Q) = 0. Note that solvability in this case implicit in Hypothesis [TJ 
Going back to the [/-equation, we find that the leading order large scale problem 
reads 

DudxxU^ + H{Uo, 0; 0) - fI{Uo)E{Uo, 0; 0) - 0, (4.10) 

where F^ is the leading order correction from the solution of the ^-equation at order 
that is V2, compare (|4.4p . The correction term F^ only depends on F^ , , and 
vanishes if either of these does, which holds in all cases analysed in the literature. 
The quadratic case. Here Vi =/= in general and so the expansion of (j4.9p at order 
contains the terms G^{Uo,Vo)V2 and duG^ {Uo,Vo)[Ui,Vi]. This means that the 
equation for Vi couples to the next orders in U and V. However, the examples of the 
quadratic case in g5]have G^{U, V) = &{U, V)V, that is, duG^{Uo, 0) = ([/q, 0) = 0. 
Hence, the terms that couple to other orders vanish and the solution can be expressed 
as in the linear case now with a correction from solving for Vi rather than V2. 

The fact that Uoj obeys a second order ODE for all j, and Ug j^rj) = uqj = wj- lies 
at roots of {F^ , G^) implies that the limiting slopes dxU^{rj±) completely determine 
the adjacent large scale solutions Unj, Uoj+i- Hence, the matching problem is local 
in the interface sequence, and the velocities cj can be viewed as a function of the left 
and right slopes (and roots of (F^,G^)) only: cj = G{dxU~{rj — ),dxU^{rj+),Uj). 
The laws of motion for first order interface dynamics in terms of G are given by p.6p . 
However, even for the simplest cases, nothing is known analytically about G 7^ 0. 
Compare ^ Nevertheless, for n = 1, the local coupling and uq = have strong 
consequences for possible interface dynamics. See fJT] 

4.3. Further comments. As mentioned in ^ the 'normal form' proposed in 
[S] for pulses and the three component model in |15j with fronts have the second 
order standard form for semi-strong interaction, and indeed the relative motion has 
velocities of order . See [51 [TBll551l34j . The abstract derivation explains why for first 
order interaction a factor UV (and not just V) is required in the nonlinear kinetics: 
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there must be must be a simultaneous root in U of the leading U- and V^-kinetics. 
Moreover, it reveals why the quadratic case UV'^ is natural for F and G. 

In case E = 0, the corrections Fl,F^ vanish and so the large scale vector field is 
the same for first and second order interaction (as in the examples from ^2.2\i . Thus 
it does not reveal the order of semi-strong interaction. 



5. Application to examples. In this section we apply the standard forms to 
immediately see how first and second order semi-strong pulse interaction arise in the 
class of models given by 



dtU = dxxU + a ~ ^iU + -fV ~ pUV^ 
dtV =e^dxxV + f3-V + UV 



2« T/ _L « _ 1/ ^ 7-n/2 (5-1) 



where U, V are scalar and a, /?, 7, p are parameters. This is a combination of the 
model from (p = ^ = ^ 0), the Schnakenberg model (/i = 7 = [TUlIll]), 
the Gray-Scott model (a = ^, 7 = /3 = [SIHTUSO]) and the Brusselator model 
(a = /i = [Hj). 

In order to sec which type of semi-strong interaction occurs, we substitute the 
scaling 

U = e''c7, V = e'V 

into (|5.ip . Upon dividing by and e'^ , respectively, this gives 

dtU =dxxU + e~'-a^^iU + je'''-V-e^'pUV^ , , 

dtV =e^dxooV + e-'(3-V + e'+'-UV^. ^ ' ' 

5.1. Second order semi-strong interaction. Comparing (|5.ip with (|3.ip . we 
infer F^ = ^ requires s — r > —1 and 2s > —1 with at least one equality, and 
r + s = — 1. Both cases yields s = r = — 1/2. which gives 

dtU ^dxxU + e-^'^a- pU + e-^-iV -e-^pilV"^ , . 

dtV =e^dxxV + e^'^P-V + UV\ ^^''^^ 

It follows that second order interaction requires a = e^l'^a. with bounded a. and 
/3 = with bounded /3 for a regular expansion, and with Vq = in mind. The 

small scale problem (|4.5p applied to this case reads 



9{Mo = 

^cgo = -ego + t'o - uo^'o- 



„ . (5.4) 



Note that the (?;, g)-equations are the same as (|2.5p . In particular only pulses exist. 

On the other hand, expanding (|5.3p in e gives Vq = and Vi = /? so that the 
large scale problem (|4.4p , which is obtained by expanding (|5.3p 1 reads 

9.xJ7o = -(a + 7/3) + M^^o, (5.5) 



in contrast to (|2.8p . 
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5.2. First order semi-strong interaction. The comparison of (j5.1[) with p.ip 
shows that first order interaction requires s — r > —2, 2s > —2 with at least one 
equality and s + r = —2. Both cases yield r = 0, s = — 1 and thus 

dtU =d,,U + a~ ^iU + e-^jV-e^^pUV^ , . 

dtV = e'^d^^V + eP -V + e-^UV'^ . ^ ' ' 

Since has a factor V this is the quadratic case from ^34.2. 21 Wc thus expand (j5.6|) 9 
in e and solve for Vi , which suggests the two solutions 



2Un \l 4t/2 Uo' 



However, the fact that J7o = at pulse locations implies that only Vi = V{' is an 
option. Substitution into the expansion of (j5.6[) yields, in contrast to (|5.5p . the large 
scale problems 



d.^Uo = -a + ^iUo - -fVf - pUo [V,- j . (5.7) 
Now we turn to the small scale problem. The application of (|4.8p to (|5.6p gives 
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^^'o = -cd^vo + vo- uiVq 



(5.8) 



Analogous to (|2.10p . for any c the invariant subspacc {vq = O^vq — 0} of (|5.8p consists 
of afhne ui{^) = ei^ + 62, 61,62 G M. Compare also (|4.8p . For the four dimensional 
flow of (|5.8p . this space is the two-dimensional center manifold of any equilibrium 
with iiQ = and has transverse eigenvalues (— c± v4 + c?) /2. Hence, it is normally 
hyperbolic with one-dimensional stable and one-dimensional unstable manifolds and 
the center manifold has the same properties as in (|2.10p . 

For 7 = 0, (0 = 1 (|5.8p is the well-known 'core problem,' first derived and nu- 
merically analysed in [30], see also [57]. In this case also ui = is invariant and ui 
convex for ui > 0. For c = 0, the existence of various even pulse solutions which 
are homoclinic in vq and asymptotically affine in ui was proven in [6]. Some first 
order pulse motion was studied in [30], and in more detail in [T^, which showed 
that asymmetric pulse-type solutions require c 7^ 0. Indeed, in contrast to (|5.4p here 
matching requires lim^_i.-|-oo '^i(C) ~ dxUo{rj±), which breaks the reversible symme- 
try if dxU[){rj—) + dxUo{rj+) 7^ 0. Therefore, motion is expected to be of order e, 
though nothing is known rigorously. 

Note that second order interaction is an asymptotic regime within first order 
interaction. Just as for the example in ^32. 21 substituting the scaling for second order 
interaction into (|5.7p and (|5.8p yields (|5.5p and (|5.4p . respectively. 



6. Lavifs of motion for second order semi-strong interaction. In this sec- 
tion we derive the leading order interaction laws in case of second order interaction 
based on the standard forms and Hypothesis [TJ To capture the leading order motion 
let us first consider the small scale PDE for v from (j3.ip for = 0, 



dfV = Dyd^^v + G*'(u, v)v + eE{u, v; e), 
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(6.1) 



We discuss the M-componcnt in i J6. II below. In the foUowing, solutions from Hypoth- 
esis [T] are considered, and for readability we suppress dependencies of u, w on j and e. 
Anticipating motion of order e^, the ansatz rj — Vjir), r = e^t is appropriate. The 
strategy is now to expand (j6.ip in e and derive solvability conditions by comparing 
terms of equal orders. This will also yield conditions that determine the shapes of 
profiles, which were already encountered in 21 

To expand the left hand side of (j6.ip . we first linearise (j6.ip in an interface pattern 
from Hypothesis [T] To leading order this yields the linear operator 

C := + G*'(wo, vo) + 9t,G*'(wo, vo)wo, 

and d^VQ G ker(£) due to translation symmetry in ^ e M. Therefore, the A^-fold 
zero eigenvalue from Hypothesis [1] stems from the relative translations of interfaces. 
Hence, substituting ^ = (x — rj)/e^ and expanding in e, the left hand side of (|6.ip 
gives 

^^v{t,C} = -ed^v^^r,+0{e^). (6.2) 

Concerning the right hand side of (|6.ip . order recovers the (wq, go)-Gquations from 
(|4.5p for c = as in p.2p . At order the right hand side of (16.11) reads 

e{Cvi+duG\ua,va)[ui,vo]+E{ua,VQ-{))) + 0{e^). 

Comparison with terms of order e from (|6.2[) implies 

-Cvi =d^vo-^rj + duG''{uo,vo)[ui,vo] + E{uo,vo;0). 

UT 

For instance in an setting there is w & ker£*, with the L^-adjoint £*. This yields 
the solvability condition 

{d^i'o -^Tj + duG''{uo,vo)[ui,vo] + E{uo,vo]0),w) = 0, 

UT 

that determines the velocity as in (|3.4|) via 

-^rj ^ -{duG^{ua,va)[ui,vo]+ E{uo,vo;0),w)/{d^va,w). (6.3) 
clr 

Here vq = vqj are the interface profiles, while uq = a-j are constant in ^, and ui = 
Uij is given in (jS.Sp . By Hypothesis [T] the translation modes are the only neutral 
eigenmodes at e = so that an A^-dimensional vector field for the second order semi 
strong interface motion is expected. However, this requires to solve for r by matching 
equations p.3p and solutions to p.2[) with solutions to (|4.4p for all j = 1,...,A'^ 
and the boundaries. We expect a globally coupled vector field, but there may be 
singularities. See, e.g., [TO] . 
On the x-scale (|6.ip reads 

dtV = e^D,d.,^V + G%U, V)V + eE{U, V- e), 

with left hand side of order . To order this recovers the condition G^{Uo, Vo)Vb = 
0. For order we expand the right hand side and obtain 

G%Uo,Vo)Vi + E{Uo,Vo;0) ^ 0, 

which yields the correction to the root Vq G {0, \/''°"'} that already appeared in ^^i^ 
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6.1. The [/-component. Recall the [/-component of p.ip for i^^ = is given 

by 

dtU = D^d^^U + H{U, V- e) + e-^F\U, V)V. (6.4) 

For second order interaction the left hand side is of order , and thus (|6.4p does 
not contribute to the leading order motion of rj directly. The arising leading order 
solvability conditions only determine the shape of the interface pattern and slightly 
extend those from With left hand side 0{e'^), expanding u,v in e on the right 

hand side yields 

Dud^^uo + e{Dud^^ui + F'{uo,vo)vq) = 0{e^). 

At order £° we find the condition d^^uo = as in (|4.5[) . and order gives 

-Du%Mi + F%uo{rj), vo)vo = 0, 

which is (|3.3p and determines ui(^) up to affine terms. 

On the x-scale, away from interfaces, Hypothesis [T] implies Vq € {0, \/f'°nt| gQ 
that at order we find 

DudM + H{Uo, Vo; 0) + F%Uo, Vo)Vi = 

as in (j4.4p . Recall that Vi was determined above. At order we abbreviate Hq = 
H{Uo, Vo;0) and F^ = F"([/o, Vq), and find the condition 

Dud^^Ui + duHoUi + d,Ho + duF^[Ui,Vi] + d.,F^[Vi, Vi] + F^iV^] = 0. 

This is closed by the condition that arises at order in the w-equation, and solvability 
is implicit in Hypothesis [TJ 

6.2. Simplifications by symmetries. All thus far analysed models where sec- 
ond order semi-strong interaction occurs have n = 1 so that C is self-adjoint. Hence, a 
natural consideration in the general case is to assume self-adjoint £. so that w = d^vo 
can be chosen. In this case we can proceed as follows for symmetric interfaces. 

6.2.1. Pulse. If the interface is a symmetric pulse then vq is even (for appropri- 
ate vq{0)) and d^vo odd. In addition, E{aj,vo;s) is even as a funtion of ^, and using 
p.3p with bj := po(0), the function po — is odd. Further, ui — wi(0)-|- Jq PoiO~^j'^C 
is even, while bjd( ~ bj^ is odd. Therefore, when including j-dependence for em- 
phasis, (|6.3p simplifies to 

-^rj = -{duG''{a.j,voj)[bj-,vo,j],d^vo^.j)2/\\d^VQj\l. (6.5) 

This is indeed a generalisation of the equations of motion for pulses in the Gray-Scott 
model |33] and the Gierer-Meinhardt model [51134]. 

The equilibrium hi = . . . = Bn = means that the pulse pattern is stationary to 
leading order if the average of the slopes bj ~ po(0) = {po{oo) + pq{—oo))/2 vanishes 
at all pulses, for instance at a symmetric configuration where all distances are equal. 
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6.2.2. Front. If the interface is a front, then the equations of motion may sim- 
pHfy as follows if there is such that vq — is odd. First, we shift v so that vq 
is odd and V^'''"' = -^^f™"* = -v,, and d^vo is even. Then and G"' have the 
factor (w — V, )(w + v*), which is even in ^ for v ~ vq. If now is odd in v then 
F^{aj,v) is odd in ^ and, by (|4.5p 9. po is even. With /?j wi(0) we then have that 
ui — Pj = Jj^^po(C)dC is odd. Including ^-dependencies, (|6.3p simplifies to 

-j^rj = -{duG'{aj,vo,j)[lij,vo,j\ + E{aj,vo^j;Q),d^vn^j)2/\\d^vo,j\\l- (6.6) 

As mentioned in fj3l in [15] = duG^ = and the fronts are odd, and indeed the 
equations of motion correspond to (|6.6p . 

Remark 1. For duG^ = and E{U,V;e) = e{U;e)V with scalar e, we have 
■^rj = since {voj,d^Voj) = 0. Hence, such E do not drive second order semi- 
strong interaction of odd fronts. Clearly, duG^ ^ is required for non-trivial (|6.5p . 

The interface coupling and velocities essentially depend on bj (averages of slopes 
of Uq at rj) and Oj (values of Uq at r^). However, the additional /3j-dependence in 
(j6.6p and that duG^ = means the right hand side of (|6.3p depends only on aj. 
Through the matching conditions, these parameters depend on the neighboring inter- 
face positions or the boundary. For the three-component model the arising globally 
coupled system of algebraic equations have been derived in |16| . For the Schnakenberg 
and Gierer-Meinhardt models these are contained in [TU] and [M], respectively, and 
can have singularities at which the manifold of pulse patterns folds. The vector field 
for interface motion in general intricately depends on the details of the model and 
boundary conditions. Numerical observations suggest gradient-like dynamics, which 
appears to hold tru^ for the three-component model of [15]. However, it seems dif- 
ficult to prove this in broader generality - for first order interaction we prove results 
in this direction below in fJT] 

The degenerate case F*' = for pulses, which was discussed after (|4.3p . implies 
that Uq is constant in time and Oj = Uo(rj), bj = dxUo{rj). Hence, the systems (|6.5p 
decouple, respectively, and each pulse moves according to the same scalar ODE to 
leading order. In particular, the reduced pulse motion is monotone, and, if global, 
each pulse converges to either an equilibrium or infinity. 

7. First order semi-strong pulse motion. In this section we study first order 
semi-strong motion of pulses for scalar large scale problems (ri = 1). We assume 
existence of smooth solutions to the reduced dynamics of p.6p . For pulse patterns 
we make the natural assumption that the states U{rj) = are all equal and thus 
may be moved to zero, as in ^32. 21 and abbreviate C(P~, P+) ~ G{P~ , P+, 0). Thus, 
pulse positions are Dirichlet boundary conditions for the (leading order) second order 
large scale problem (|4.10p . whose solution inbetween pulses (if it exists) is therefore 
generically determined by the pulse positions alone. 

According to p.6p the reduced first order semi-strong dynamics on the time scale 
T ^ et with interfaces at ri < . . . < rjv, boundary conditions at tq, r^^i € RU{±c»}, 
and J{U) H{U, 0; 0) - FI{U)E{U, 0; 0) from (liH))) then reads 

DudxxU = -J{U) , X eUj=o^Nirj,rj+i) 

Uirj) - (7.1) 
^ C{dxU{r,-),dxU{r,+)), 



P. van Heijster. Personal communication. 
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where j ~ 1,...,N. For bounded D we assume separated or periodic boundary 
conditions, and in unbounded directions convergence to constant states. Note that U 
here is Uq in the notation of the previous sections. In this section subindices of U, P 
are not related to the expansion in s. 

The function C{P^ ,P^) is an essential part of the dynamics and effectively pa- 
rametcrises the manifold of pulse patterns. As mentioned, even for basic examples 
nothing is known analytically about interfaces for C 7^ 0. For P~ + P"*" = and 
C = 0, existence results are given in [5], but these do not cover the numerically 
observed folding of C as P+ increases. See [l0l[26l[30] . 

We will show that under suitable assumptions on C and J the first order iV-pulse 
motion is gradient-like with respect to various geometrically meaningful Lyapunov- 
functionals: the largest interpulse amplitude and distance decay in time, while the 
smallest of these increase. This severely constrains the leading order pulse dynamics. 

In preparation, we define the following key quantites. 

Definition 7.1. Let U{x) solve DuUxx + J{U) = with initial conditions U{0) = 
0, and dxU{0) = P. Let A±(P) =/= be the smallest positive, respectively largest 
negative value of x such that U{x) ~ 0, and set A±(P) := ±00 correspondingly if 
there is no such point. 

Note that due to (j7.1|l 9. if C/ 7^ in (r^ , rj^i) then 

r,+i - rj - A+idxUir,+)) = A_ (9,C/(7 ^-+1-)). 

Definition 7.2. Let U{x) solve DuUxx + J{U) = with initial conditions ?7(0) = 
U, and dxU{0) ~ 0. Let A*{U) ^ be the smallest positive value of x such that 
dxU{x) ~ 0, and set A*{U) := 00 if there is no such point. 

We shall show that the aforementioned decay and growth of the largest and small- 
est amplitudes between pulses occurs, as long as the pulse positions lie in the following 
set. For the distances a stronger assumption is required; sec Corollary 17.61 Fig- 
ures [7?ll ESI illustrate the results. 

Definition 7.3. For n = 1 let Vn C M'^ be the set of pulse positions r{T), 
ri(r) < ... < rN^r), which solve (j7.ip for t in an open set with hounded U{t,x), 
x € D, that is continuously differ entiable in t and such that the following hold. 

1. U{t, a;) > for x 7^ rj (r), 

2. A*(C/) is bounded and grows strictly for U in a neighborhood ofmax{U{T,x) : 
rj (t) < X < rj^i (r)} for all j = I, . . . , N — 1, and for j = and j ~ N if ro 
and rjsr+i is bounded, respectively, 

3. sgn{CiPr,P+)) = sgn(P,r + P+). 

A priori the set Vn could be empty, but based on Remark [2] below and numerical 
evidence in TO'j'W this is for instance not the case for the relevant solutions of the 
'core problem' (15. 8p . However, there C(P~,P"'") need not be well-defined globally in 
time as the dynamics might be driven over the aforementioned fold of the graph of C. 
Another possibility for the PDE dynamics to leave the slow manifold of pulse patterns 
is a transverse instability of the manifold, beyond which (|7.ip is not meaningful. 
Numerical evidence in [10] shows that this can indeed occur. 

Remark 2. The monotonicity requirement of Definition ] 7. ,'J]f 2 ) holds for solutions 
U of dLl]) if J'{U) < and J{U)U > for x G (r^-, rj+i). (The proof is given at the 
end of this section.) This is always the case near saddle points (which means relatively 
large pulse distance). The known concrete models from ^ for first order semi-strong 
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min A^(r) 




max Al (r) 

Fig. 7.1. Illustration of the Lyapunov Junctionals and interface dynamics for periodic D and 
N = 4. In this case there are no further local extrema of M{U), Ai{r). 

interaction have J of the form J{U) = ei — e2U so that the monotonicity holds for 
all relevant solutions i/ 61,62 > 0, and U < 61/62. 

The general monotonicity problem of A* is closely related to the monotonicity 
of the period function which has been extensively studied in the literature. See, 
e.g., [13j and the references therein for a recent account. A fairly practical criterion 
for monotonicity given in [5] is that V/J"^ be convex. 

Remark 3. The sign condition U > except at pulse positions in is not 
necessary, but holds for the model class in J3 The case U < can be treated analo- 
gously. In fact, if U lies on a periodic orbit, then an even number of sign changes of U 
between pulses can be removed by shortening D without changing the local dynamics. 

We define the set of local maxima of a large scale solution by 
Af{U) := {U{x) : U (x) is a local maximum for x £ D, x ^ rj, j — 1, . . . , N} , 
and with r = (ri, . . . , r^r) the distances between pulses by 

M{r) :={d, =r,+i-r, :j=0,...,iV}. 

For X € (ro, ri) U (r^r, rjy^i), a technical issue is to identify an adjusted maximum 
and distance that are suitably comparable to those inbetween pulses. Therefore, we 
first consider periodic D so that j S Z mod + 1. 

Definition 7.4. Let j ^ 0, . . . , N. 

1. Let U* := sup{C/(x) : x € We call U* a local maximum of N{U) 
ifU* = U* for some j and U* > U*_^, U* > U*^^ or U* > U*_^, U* > U*^^. 

2. We call d* a local maximum of M{r) if d* = dj < 00 for some j and if 
d* > dj-i, d* > dj+i or d* > dj-i, d* > rf^+i. 

3. Local minima are defined analogously in each case. 

The main results are the following. 

Theorem 7.5. Assume n = I, periodic D and r = r(r) e 2?jv with associated 
U = U{t, x) . For each j ~ 0, . . . , N the following holds. 

1. IfU* is a local maximum of M{U) then -j^U* < and -^dj < 0. 

2. IfU* is a local minimum of M(U) then j^U* > and -^dj > 0. 
3- (^) 4^U*^0^ (n) Arf^. = ^ i^n) uf^ U*_, = U*+l 

Corollary 7.6. Under the assumptions of Theorem \7.5\ it holds that maxA/'(C/) 
and — minA/'(?7) are strict Lyapunov functionals in the sense that these are either 
constant with constant r, or the functionals decay strictly in t until r lies in the 
boundary ofT>N (which may never happen). If in addition A* is bounded and strictly 
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U 
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Fig. 7.2. Sketch of the phase space geometry. Thin lines are trajectories (separatrices in 
grey), bullets equilibria. The pictures extend by reflection about the U-axis and reversing arrows on 
trajectories, (a) For the proof of Theorem \ 7. 5\ The thick gray arrows indicate the dynamics with 
respect to t for a local maximum Ut G Af{U). (b) For the proof of Remark\^ 



monotone increasing on the interval (minA/'([/), maxA/'(C/)), then also maxA^(r) and 
— minA^(r) are Lyapunov Junctionals in this sense. 

Proof. (Theorem [7j5]). For n = 1 the scalar large scale problem is reversibly 
symmetric (from reflection symmetry in x) with phase space (U,P), P = dxU. For 
the given solution we denote {Uj,Pj) := {U,dxU) on [rj,rj+i]. 

Since Uj > in (rj,rj+i) and Uj{rj) = Uj(rjj^\) — 0, the reversible symmetry 
about the t/-axis implies that a unique intersection of {Uj,Pj) with the [/-axis occurs 
at U* e M{U). In particular, dj = A+(Pj(rj+)) and Pj{rj +dj/2) = 0, which means 
U* = U,{r,+d,/2). 

Such solutions are ordered in the following sense. Consider solutions [/ > to 
(HUli with t/(0) = 0, dxU{Q) = P* and A+(P*) < oo. Since the phase space is two- 
dimensional, the maxima U* ~ [/(A+(P*)/2) are strictly increasing in P* as long as 
A+(P*) < oo. See Figure O 

For U* a local maximum of A/'(C/), set P^ := P{rj±). The ordering for interpulse 
profiles implies that -Pj' < P+ and -Pj'+i > Pj+i, or -Pj^ < Pj^ and -Pj'+i > 
Pj+i- By Definition I7.3f 3') we have rj > and rj+i < or rj > and rj+i < 0. 
Therefore, 

^(^.-^.+i)<0, (7.2) 

and monotonicity of A_|- in Definition 17. 2) implies that j^U* < 0. Analogously we 
find that -^U* > in case of a local minimum, which proves the first two parts of 
the theorem. 

For the third part, it follows from the above that -^U* ^ implies U* ^ 
for k = j — 1 or k = j + 1. Since U* ~ implies that Uj and Uk lie on the same 
trajectory, equivalence of (i) and (iii) follows from Definition 17.3( 3). Equivalence of 
(i) and (ii) is a direct consequence of Definition 17.3( 21. □ 

Remark 4. Due to the global ordering of trajectories noted in the proof, there are 
many more such Lyapunov-functionals: in (U, P)-.space, the intersection points with 
any fixed line through the origin in the positive quadrant move towards the origin. See 
Figure\77S\ 

Proof. fCorollarv l7.6[) . The first part immediately follows from Theorem 1 7. 5 1 
Concerning the second part, it follows from the global monotonicity assumption 

that if dj is a local extremum of Ai{r), then U* is a local extremum of Af{U). Now 

application of the proof of Theorem 17.51 implies the claim. □ 
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Fig. 7.3. Illustration of the Lyapunov functions and interface dynamics for bounded D with 
separated boundary conditions and N = 5. Dashed lines mark the extensions beyond D. The 
extension on the left generates local minima of M{U) and M{r), and that on the right local maxima. 



Proof. (Remark 12). Generally, [/(r^) = U{rj+i) = and sgn{U{x)) = ±1 in 
(rj , Tj+i) requires sgn(J(?7)) = ±1 at the critical point U* . Hence, J{U) ^ requires 
sgn(J(f/)) = sgn([/) for solvability of (j7.1|) . We consider U,J{U) > 0; the negative 
case follows by a symmetric argument. Hence, wc have dxP = —J{U) < so that P 
decays strictly. 

Let (C/i,Pi), (C/2,P2) be two solutions to DuU^:, + J{U) = with < Pi(0) < 
P2(0), Ui{0) = C/2(0) = 0. Let > be smallest so that Pj{xj) = 0. We need to 
show that xi < X2. 

The strict decay of P2 implies that there is smallest > so that P2(a;*) = Pi(0). 
See Figure rr^ b). From the ordering of trajectories noted in the proof of Theorcm l7.5l 
we have that Pi(?/i) = P2(?/2), imphcs Ux{yi) < U2{y2) for yj e (0,a;j]. Now J > 
and J' < gives > -J(C/2(y2)) > -j\Ui{yi)). Hence, d^Pj = -J{Uj) implies 
Pi{x) < P2{x* + x) for all x G [0,a;i] and so 2:2 > xi. □ 

7.1. Separated boundary conditions and unbounded D. 

To make the boundary segments Uq and Un comparable with the intcrpulse seg- 
ments, we define := ri + A^{dxU{ri~)) and rjv+i := fjv + A+{dxU{r]y+)). We 
further adjust Definition 17.41 and the definitions of Af{U), Ai{r) to non-periodic do- 
mains as follows: (1) replace ro with fg and rpf+i by rjv+i, (2) set U^i = Uq, d-i = do 
and ^ ^iv+i' '^JV+2 = djy+i. 

Theorem 7.7. Assume non-periodic D. Theorem \ 7. 5\ holds with respect to the 
adjusted M {U) , AA{r). Moreover, if {rQ^rN+i) is hounded, then Corollary \ 7. 6'| holds 
for the Lyapunov functional^ defined on [^Oi^A'+i] with the adjusted M [U) , A^(r). 

Proof. With the adjustments the proofs are the same as those of Theorem 17.51 
and Corollarv 17.61 respectively. □ 

Theorem 7.8. Assume unbounded D, r e 2?jv and that U converges to constant 
states in unbounded directions. Let j £ {0, + 1}. // \rj\ = 00 then sgn(^rj) is 
constant, and nonzero if N > \. If in addition J has a unique saddle in {U > 0} and 
N >l then sgn(^rj) = sgn(rj). 

Note that the models in fJ5]have a unique saddle if > 0. 

Proof. We consider only ro = —00 as ttv+i = 00 follows by symmetry. The 
convergence to an equilibrium implies that this must be a saddle point (since n = 
1) and Uo{x) := U{x) for x < rg equals its unstable manifold up to the point 
{Uvi{ri),dxU{){ri—)). Therefore, dxUQ{ri — ) is constant in time. The stable mani- 
fold is the reflection of the unstable one about the [/-axis, and forms a separatrix 
in the planar phase space. Hence, for all r either dxUQ{ri+) 7^ —dxUo{ri — ), or 
dxU{){ri+) = —dxUo{ri — ). The latter implies that Ui lies in the stable manifold of 
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the saddle, which requires = 1 and = 0. 

Concerning the former, following the arguments of the proof of Theorem 17.51 
(trajectory ordering and the definition of Pat) shows that sgn(^ri) is constant in 
time. See Figure [721 a). 

For > 1, the point {U{ri), dxU{ri+)) cannot lie in the stable manifold of 
the saddle point, and thus dxUo{ri—) =^ —dxUo{ri+), which implies j^ri 7^ by 
Definition ESIJS). 

If the saddle is unique in {[/ > 0} then [/(ri) = together with ?7 > for x < ri 
implies that the saddle has tZ-coordinate [/* > 0. Uniqueness of the saddle and that 
its stable and unstable manifolds form separatrices imply the following: trajectories 
with U > which intersect {U = 0} intersect the J7-axis in the interval [0, U^,). See 
Figure Ella). As in the proof of Theorem O it follows that < if iV > 1. □ 
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